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Cd The polymerization of actin via branching at a ceU membrane containing nucleation-promoting 

^^ factors (NPFs) is simulated using a stochastic-growth methodology. The polymerized-actin dis- 

OO 

CN tribution displays three types of behavior: a) traveling waves, b) moving patches, and c) random 

r~| fluctuations. Increasing actin concentration causes a transition from patches to waves. The waves 

I and patches move by a treadmilling mechanism which does not require myosin II. The effects of 
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downregulation of key proteins on actin wave behavior are evaluated. 



The dynamic behavior of the intracellular protein actin in cells is crucial in directing cell 
shape changes and migration. Actin occurs in both monomeric ("G-actin"), and polymerized 
( "F-actin" ) forms. F-actin consists of polar semiflexible filaments which grow preferentially 
at their "barbed" ends. The filaments assemble into supramolecular structures which provide 
forces to move cells or form dynamic protrusions along the cell periphery. A widely observed 
type of supramolecular structure is a branched dendritic network, where new filaments grow 
as branches on existing ones [H |2]. Several investigations have shown that F-actin in cells 
displays spontaneous dynamic behavior including traveling waves and patches pHH]. The 
waves and patches in Dictyostelium typically move at about 0.2 ^m/s. A closely related 
phenomenon is Hem-1 waves in neutrophils, where Hem-1 acts upstream of actin polymer- 
ization [7]. Actin waves, when they impinge the membrane, cause cellular protrusions to 
form ISJ. Therefore they may be involved in oscillatory and/or random motions of the cell 
membrane, which allow a cell to explore its environment. Thus understanding the origins 
of actin waves can help clarify the mechanisms by which cells determine key properties of 
their motion. 

The molecular mechanisms underlying actin patches and waves are not well known. How- 
ever, a partial molecular inventory of the waves has been obtained [S], including Arp2/3 
complex, which nucleates new filaments by branching, myosin-lB, which may link the actin 
network to the membrane, and coronin, which disassembles actin networks. The presence 
of Arp2/3 complex suggests that patches could consist of dendritic actin networks; patches 
in yeast do in fact consist of such networks |9] . 

Actin waves and patches might result [3l [T014T2] from a react ion- diffusion mechanism 
|13j . based on spatially varying concentrations of an "activator" and an "inhibitor". The 
activator grows at the front of the wave by positive feedback, and spreading (for example by 
diffusion) drives the wave forward. The inhibitor then builds up, suppressing the activator 
at the rear of the wave. Actin waves and patches could contain positive feedback either in 
F-actin ("actin- first"), in membrane-bound filament nucleation-promoting factors which act 
upstream of actin polymerization ("NPF-first"), or both. Ref. [11] treated a coarse-grained 
model of F-actin dynamics with delayed membrane-induced inhibition and positive-feedback 
of F-actin on itself, which could implicitly contain contributions from NPFs. Assuming 
spontaneous polarization of actin filaments and diffusion-like spreading of F-actin led to the 
formation of patches and waves in appropriate parameter ranges. Refs. [3| and [12] found 



spontaneous waves of F-actin and the NPF Hem-l in models based on NPF cooperativity 
and inhibition of Heni-1 by F-actin; Ref. [12] also found actin patches under some conditions. 

This Letter uses a detailed 3D dendritic network model to treat the spatial dynamics of 
F-actin, coupled to a NPF which shuttles between the cytoplasm and the membrane. This 
allows a treatment of actin waves and patches based on well-established molecular mecha- 
nisms, permits exploration of the possible dynamic phases of F-actin, and enables prediction 
of the dependence of F-actin dynamics on key biophysical rates and protein concentrations. 
The main finding is that known molecular features of dendritic actin nucleation embody 
mechanisms of positive feedback, spreading, and delayed negative feedback which generate 
traveling waves and patches. This model differs from that of Ref. [11], in that polarization 
of filament orientations is not required. Three broad types of F-actin dynamics are found: 
1) traveling waves, 2) moving patches, and 3) random fluctuations with occasional moving 
patches. 

The simulation method extends one used previously to treat dendritic actin networks 
[m |T5]. I treat network growth on a 3fim x 3fim piece of membrane. Explicit coordinates 
are stored for each filament. Polymerization, filament nucleation by branching, and capping, 
which prevents barbed-end polymerization, are treated stochastically. New filaments are 
generated near the membrane as branches on the sides of existing filaments, and by a slower 
non-branching nucleation process. Branching is allowed only within 10 nm of the membrane, 
so that most branches form near filament tips. The method is extended by including the 
following effects (see Supplementary Material for more details): 

(i) Filament severing, which is the primary actin disassembly mechanism in the model. 
It is also an important actin filament nucleation mechanism [TU] . 

(ii) Attractive filament-membrane forces, with binding energy Ef,, required because the 
actin network is attached to the membrane [8J. 

(in) Translational and rotational motion of filaments, which are important for filament 
detachment from the membrane and for "docking" of filament clusters to the membrane. 
Because this is very computationally demanding, I simplify the calculation by treating con- 
nected dendritic filament clusters as rigidly moving units. The translational cluster motion 
has a deterministic part from membrane forces, and a random part from Brownian motion. 

(iv) Membrane attachment, detachment, and diffusion of NPFs. Detachment reduces or 
eliminates NPF activity, and F-actin enhances the detachment rate [7j. Taking the mem- 



brane as the x-y plane, the attached (active) NPF density na{x,y) and spatially averaged 
detached activator density na satisfy 

(JT) 

-7^ = -kdetF{x, y)na + KttUd + -DmcmbV^na 

at 

dt ~ ^ dt ^' ^ ^ 

where fcdet, ^att are constants, -Dmemb is the attached NPF diffusion coefficient, () denotes 
2D spatial averaging, and F{x,y) is the 2D density of F-actin (see Supplemental Material). 
Detached-NPF diffusion is assumed fast enough that nd{x,y) is constant. Because dimer- 
ization greatly accelerates the activity of attached NPFs [17J, I take the branching rate for 
a filament impinging the membrane at {x,y) to be proportional to na{x,yY. 

The parameters (see Supplemental Material) are taken from experiments where available 
and otherwise are chosen to obtain a realistic network structure. The on- and off-rates, the 
severing rate, the branch dissociation rate, and -Dmemb are obtained from experimental data. 
The value of [A] is varied within biologically reasonable limits, and the branching rate and 
[CP] are chosen to give reasonable network structures; i^b, ^att, ^det, and the non-branching 
nucleation rate fcnuc are treated as independent variables. 

The simulations begin with only G-actin present and are run out to 400s. Polymerization 
peaks and later reaches a steady-state dynamic regime, which is my focus. Fig. 1 shows 
snapshots of the 3D actin network and the simulated fluorescence intensity (with a spread 
of 200 nm), under conditions favoring wave formation. I define a wave to be a moving 
feature having a large aspect ratio, while features with smaller aspect ratios are patches; 
a quantitative description is given in the Supplementary Material. A wave of polymerized 
actin moves across the membrane at a speed of about 0.2fi'm/s. The wave's shape changes as 
it moves, due to the randomness of polymerization and detachment. Although the F-actin 
features are mainly wavelike, switches to a patchy state and back are sometimes seen in the 
simulations (see end of Video 1). 

At lower [A], a patch initiates (left frame. Fig. 2) as a localized focus of actin polymer- 
ization. It remains stationary for about 4s and then moves to the right at about 0.15 fim/s, 
maintaining a roughly constant size. The patch forms because G-actin becomes depleted 
before a complete wave can be formed. This result is consistent with previous studies of 
reaction-diffusion systems \TS], which showed that global inhibition such as that mediated 
by G-actin depletion causes transitions from waves to patches. The reduced [CP] pictures in 




FIG. 1: 3D network structure (left) and simulated fluorescence (right) for actin wave (red) propa- 
gating across three- micron square area of membrane (green) . Time difference between frames is 6 
sec. Parameters are: k^tt = 0.025s-^ k^et = 0.009^J.M-^s-^, E^, = 2.8/cbT, knuc = 0.003/iM-ls-^ 
[A\ = 22.5/iM, and [CP] = 0.25^M. See Supplemental Video 1 for complete time evolution. 



Fig. 2 show a more random distribution of F-actin, with some patchy features which move 
and give the appearance of broken wave fronts. 

The motion of the waves and patches results from three factors operating jointly: 
Positive feedback of F-actin. This is casued by autocatalytic branching from existing 
filaments, which by itself gives linear positive feedback [19] that is confirmed by in vitro 
polymerization dynamics [20j. Branching is modulated by filament detachment from the 
membrane. If E\, is small, most single filaments leave the membrane before they branch, so 
that formation of a dendritic cluster consisting of several filaments, by a statistical fluctua- 
tion, is required before continued growth occurs. The requirement for such a critical cluster 
at small E\^ values results in a sudden onset of polymerization, which favors wave/patch 
formation. This effect is similar to that of nonlinearity in the positive feedback. 

Spreading of F-actin. Fig. 3 shows a filament growing from a patch by polymerization 
and nucleating several branches, thus providing new polymerization nucleui and spreading 
the patch to the right. This mechanism does not require the polarization of filament ori- 




FIG. 2: Simulated fluorescence of patches obtained with parameters of Fig. 1 but with [A] reduced 
to IbfiM (top row) or [CP] decreased to 0.10/iM (bottom row). Images are 12 sec (top row) and 
18 sec (bottom row) apart. See Supplemental Videos 2 and 3 for complete time evolution. 



entations suggested in Ref. |TTj, since a broad range of filament orientations, including the 
vertical orientation, can generate new branches in the direction of motion (see Supplemen- 
tary Material). However, the patches and waves are polarized in that there are more free 
barbed ends at the leading edge. 

Delayed negative feedback. This effect is, as mentioned above, grounded in exper- 
imental observations. The delay has a large effect on the F-actin dynamics when 1/fcatt 
exceeds the lifetime of polymerized actin, about 5 sec in the present simulations. 

These three factors can give rise to traveling waves and, when combined with global 
inhibition such as that arising from G-actin depletion, patches [18]. Small values of the 
non-branching nucleation rate k^nc are also required in order to avoid "parasitic" nucleation 
interfering with wave or patch formation. 

This model of actin waves and patches makes several experimentally testable predictions: 
1) Increasing G-actin concentration leads to a transition from stationary patches, to moving 
patches, to waves (cf. Figs. 1 and 2). This prediction is consistent with experiments [6] 
that treated Dictyostelium cells with latrunculin, which depolymerized actin by sequestering 




FIG. 3: Closeup of edge of actin patch at top of Fig. 2 moving to right. Frames are 0.2 sec apart. 



G-actin monomers. The latrunculin was then removed from the growth medium. During 

the resulting recovery of polymerizable G-actin, the F-actin distributions showed transitions 

from stationary patches, to moving patches, to waves. The patch-to-wave transition may 

well result from changes in [A], as in the present model. Because intracellular latrunculin 

is strongly bound to G-actin, the latrunculin exit required for G-actin recovery could take 

minutes or longer. The coexistence of waves and patches seen in Ref. [6] is also consistent 

with the wave-patch switching seen in the simulations. 

2) Colliding waves annihilate each other, since each depletes the NPF required for the other 

to persist. This is consistent with the experiments of Ref. [8] 

3) The waves and patches move by a treadmilling mechanism independent of myosin II (cf. 

Fig. 3). This is consistent with photobleaching experiments |8], in which a bleached spot 

in the F-actin remained stationary during wave motion and myosin II was not required for 

wave motion. 

4) Arp2/3 complex is uniformly distributed throughout the wave (see Supplemental Mate- 
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FIG. 4: Total number of polymerized subunits A'^p vs time: a) Baseline parameters as in Fig. 1; b) 
[CP] decreased to O.lO/uM.; c) knuc increased to 0.03fiM^^s^^; d) E'b increased to A.Bk-BT; e) /catt 
increased to 0.05s~^ 

rial), consistent with the fluorescence measurements of Ref. |8]. 

5) Uncapped barbed ends are polarized toward the front edge of patches and waves. This 
prediction could be tested by standard techniques for measuring the spatial distribution of 
free barbed ends. 

6) The actin filament orientations are isotropic rather than being biased along the direction 
of motion. This prediction could be tested by cryo-electron microscopy measurements of 
wave or patch structure. 

7) Downregulation of capping protein, increased nucleation activity in the membrane (as 
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from forinins), increased filament-nienibrane attachment strength (controlled by the protein 
VASP [23]), or more rapid reattachment of NPFs to the membrane, will prevent wave 
or patch formation. Fig. 4 shows the effects of such interventions on the dynamics of 
the number of polymerized subunits Np, used here as an approximate descriptor. Large 
fluctuations in A^p (as in frame a) correspond to traveling-wave or moving-patch states 
seen in the simulation videos, while smaller fluctuations correspond to randomly fluctuating 
states. All of the interventions reduce wave behavior. A recent study [21] suggested that 
CP enhances branching in vitro, so that reduced CP levels might reduce polymerization. 
However, in Dictyostelium reduced CP enhances polymerization |22], consistent with the 
assumptions made here. If CP enhances branching, the effects seen in Fig. 4b will be 
reduced. 

The main conclusion of this work is that essential mechanisms needed for formation of 
traveling waves and patches of F-actin follow from known properties of actin filament nucle- 
ation at membranes. Thus the predictions are not limited to the particular mathematical 
model employed, but rather follow from the general mechanisms of F-actin positive feedback, 
F-actin spreading, and delayed negative feedback. The first two of these follow immediately 
from the nature of branching nucleation at membranes, while delayed negative feedback 
follows from plausible interactions mechanisms between F-actin and NPFs. In combination, 
these mechanisms lead to dynamic features closely paralleling those seen in experiments. 

I appreciate informative discussions with Alexander Mikhailov, Martin Falcke, and 
Karsten Kruse. This work was supported by the National Institutes of Health under Grant 
ROl GM086882. 
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